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Different theoretical approaches for the thermodynamic properties and the equation of state for 
multicomponent mixtures of nonadditive hard spheres in d dimensions are presented in a unified way. 
These include the theory by Hamad, our previous formulation, the original MIX1 theory, a recently 
proposed modified MIX1 theory, as well as a nonlinear extension of the MIX1 theory proposed in 
this paper. Explicit expressions for the compressibility factor, Helmholtz free energy, and second, 
third, and fourth virial coefficients are provided. A comparison is carried out with recent Monte 
Carlo data for the virial coefficients of asymmetric mixtures and with available simulation data for 
the compressibility factor, the critical consolute point, and the liquid-liquid coexistence curves. The 
merits and limitations of each theory are pointed out. 



I. INTRODUCTION 

Nonadditive hard spheres represent a versatile model 
to study various real physical systems. These include al- 
loys, aqueous electrolyte solutions, molten salts, rare gas 
mixtures, and colloids. In these systems homocoordina- 
tion and heterocoordination may be interpreted in terms 
of excluded volume effects due to nonadditivity of the 
repulsive (hard-core) part of the intermolecular poten- 
tial and so, for instance, the occurrence of liquid-liquid 
demixing in real systems may be linked to a binary hard- 
sphere mixture with positive nonadditivity, while nega- 
tive nonadditivity may be invoked to explain chemical 
short-range order in amorphous and liquid binary mix- 
tures with preferred heterocoordination. On the theoret- 
ical side, prototype models of nonadditive hard-sphere 
mixtures such as the Widom-Rowlinson model^ or the 
Asakura-Oosawa model^ have been very useful to gain 
insight into interesting physical aspects such as fluid-fluid 
phase transitions and the nature of depletion forces. 

A few years ago, in a paper— where a rather thorough 
review of the theoretical and simulation work on non- 
additive hard-sphere mixtures was provided, we intro- 
duced an equation of state of multicomponent nonaddi- 
tive hard-sphere mixtures in d dimensions. Such an equa- 
tion of state results from a natural extension of the one 
we had earlier proposed for additive hard spheres^ has 
an explicit (simple) density dependence, and by construc- 
tion leads to the exact second and third virial coefficients. 
In the case of d = 3, in the same paper we compared the 
predictions for the compressibility factor corresponding 
to our proposal with those of the proposal by Hamad, 5-8 
which shares some characteristics with ours, and avail- 
able simulation results for various binary mixtures^— 
We also compared the predictions of the fourth and fifth 
virial coefficients arising from the above two theoretical 
proposals and the then available Monte Carlo resultsi 14 i 15 
The restriction in the comparison only to Hamad's ap- 
proach was justified then by the fact that Hamad had 
already proved that his proposal was superior to other 
theories, including the so-called MIXI theory originally 



due to Melnick and SawfordJ^ 

Recently, Pellicane et alJ^- have reported new evalu- 
ations of the fourth virial coefficient of a binary nonad- 
ditive hard-sphere mixture covering a wide range of size 
ratios and values of the nonadditivity parameter. Also 
recently, Paricaud 18 has proposed a new equation of state 
for nonadditive hard-sphere mixtures which is based on 
and corrects one of the deficiencies of the MIXI theory, 
namely the fact that MIXI does not lead to the correct 
second virial coefficient. These two recent papers serve 
as a motivation for the present contribution. On the 
one hand, we want to see to what extent the conclusions 
drawn from the analysis carried out in Ref. 3; are still 
valid in view of the new available data. On the other 
hand, we will also introduce a (new) nonlinear extension 
of the MIXI theory. As an extra bonus, we will write all 
the theoretical expressions in a unified language which 
will hopefully make the comparison much easier. 

The paper is organized as follows. In order to make it 
self-contained, in the next section we provide the neces- 
sary background for the later development. Section IIIII 
provides the explicit expressions for the contact values of 
the radial distribution functions, compressibility factors, 
Helmholtz free energies, and second, third, and fourth 
virial coefficients as given by the original MIXI theory, 
Paricaud's modified MIXI theory, Hamad's theory, and 
our earlier proposal. A nonlinear extension of the MIXI 
theory is also introduced here. In Sec. lIVI we compare the 
numerical values of the composition-independent virial 
coefficients, compressibility factors, and liquid-liquid co- 
existence curves for a variety of cases with available 
Monte Carlo data. The paper is closed in Sec. |V] with 
some concluding remarks. 

II. GENERAL BACKGROUND 

We consider an A^-component mixture of nonaddi- 
tive hard spheres in d dimensions. Let ov,- denote the 
hard-core distance of the interaction between a sphere 
of species i and a sphere of species j. If the diame- 
ter of a sphere of species i is <ii = an, then ov, = 



2 



\{(Ti + Cj)(l + Ajj), where Ay > —1 is a symmetric 
matrix with zero diagonal elements (A ri = 0) that char- 
acterizes the degree of nonadditivity of the interactions. 
In the case of a binary mixture (N = 2), the only non- 
additivity parameter is A = A12 = A21. 

The compressibility factor Z = p/pksT of the non- 
additive mixture, where p is the total number density, 
p is the pressure, T is the temperature, and ks is the 
Boltzmann constant, is given by 

N 

Z(p, {x k }, {vice}) = 1 + 2 d ~ 1 v d p XjXjafj 

X9ij(p,{ x k},Wke}), (2.1) 

where v d = (it / A) d / 2 /T(l + d/2) is the volume of a d- 
dimensional sphere of unit diameter, Xi — pt/p is the 
mole fraction of species i, pi being the partial number 
density of particles of species i, and g%j{p, {x^}, {ate}) = 
gij (p) stands for the radial distribution functions at con- 
tact. Unfortunately, no general expression is known for 
gij(p), but it may formally be expanded in a power series 
in density as 

N N 

9ij(p) = 1 + 

%k%£Ckl;ij 
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-o(p 3 ), 



(2.2) 



where the coefficients c^ij, Ckt-,ij, ■ ■ ■ are independent of 
the mole fractions but in general depend in a non trivial 
way on the set of diameters {o^}- To our knowledge, 
only the coefficients linear in p (i.e., c^ij) are known 
analytically for d < 3. This formal series expansion in 
the number density, Eq. (|2.2p . when substituted into Eq. 
(|2.ip . yields the virial expansion of Z which we write in 
the form 



Z(p) = l + Y d p n B n+l 
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= 1 + p £ XiXjBij + p 2 ^2 XiXjXkBijk 

i,J=l i,j,k—l 
N 

+p 3 ^2 XiXjXkXeBijke + 0(p 4 ). (2.3) 

i,j,k,&=i 

Here B n is the usual nth virial coefficient of the multi- 
component mixture, which is a polynomial of degree n in 
the mole fraction, Bij... being composition-independent 
coefficients. In terms of the coefficients Ck-ij and Ckt-ij, 
the composition-independent second, third, and fourth 
virial coefficients are given by 



Bij = 2 d ^v d at, 



(2.4) 
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Along the path we have taken, the different theories 
for mixtures of nonadditive hard spheres in d dimensions 
may be related to different proposals for gij(p). In the 
next section we provide the explicit expressions for the 
approximate proposals that we will consider in this pa- 
per, including a new nonlinear extension of the MIX1 
theory. 



III. SOME APPROXIMATE THEORETICAL 
DEVELOPMENTS 

A. MIX1 approximation 

The original MIX1 approximation^ which we will in- 
dicate with a superscript M, is equivalent to 



d M/ \ 



(3.1) 



where gf dd (p) are the contact values of the additive mix- 
ture and 



F M = dA- 



(3.2) 



Inserting Eq. (I3.ip into Eq. (|2.1j) one gets 
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Z M (p) = Z* dd (p) + b 2 v dP J2^ 



'1 < u j 



d 



(3.3) 



with Z add (p) the compressibility factor of the additive 
mixture with the same sets of mole fractions {xk} and 
diameters {uk}- The Helmholtz free energy per particle 
in the MIX1 theory is then 

aM (p) , V- , / x<A «cx d (p) . 

' - -1 + 22 Xi In (x iP \f) + e * ^' + b 2 v dP 
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Y<?9t? d (p), (3.4) 



where Xi is the de Broglie wavelength of particles of 
species i, a^ dd (p) is the excess Helmholtz free energy per 
particle of the additive mixture and, for convenience, we 
have identified 2 d ~ 1 with the reduced second virial coef- 
ficient in the one-component d-dimensional hard-sphere 
fluid 62- The second, third, and fourth virial coefficients 
of the mixture are in turn given by 
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(3.7) 



In Eqs. (JSU) and (|37f) . and cj$^ correspond to the 
coefficients in the expansion of gf dd (p) in powers of the 
number density. Note that the second virial coefficient 
of the mixture in the MIX1 theory is not exact [compare 
Eqs. (|2.4p and (|3.5[) ]. except to first order in Ay. This 
problem can be traced back to the fact that, according 
to Eq. rf3~Tjl . 



1 + Y M 
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or, equivalently, 
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In this way, instead of Eq. ([37 



(3.11) 

we have 

linip^o 9fj M {p) — 1 and thus the second virial coefficient 
becomes exact. Otherwise, the third and higher virial 
coefficients are still approximate. In particular, the third 
and fourth virial coefficients are given by Eqs. (|3.6[) and 
(1377)1 . respectively, with Y?f ->• Y£ M . 



C. Hamad's proposal 

Hamad's approximation^— denoted by a superscript 
H, consists of proposing the following ansatz 



gfM = 9 purc (y)\ y=v x» > 



(3.12) 



where ff purc (y) is the contact value of the radial dis- 
tribution function of the one-component (/-dimensional 
hard-sphere fluid at the packing fraction y, r\ = Vdp{c d ) 
is the packing fraction of the mixture (with (a m ) — 

Y^i=i x i a T)i an d Xjj wm be specified later. From Eq. 
Q3.12[) it follows that the virial expansion of gy (p) is given 
by 



9?i(p) = l + 



K+2 



n=l 



'v dP {v d )x*y 



(3.13) 



where b n is the reduced nth virial coefficient of the one- 
component d-dimensional hard-sphere fluid. In particu- 
lar, comparing Eq. (|3. 13|) with Eq. (|2.2I) . one gets 
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,,H 



(3.14) 



This is remedied by Paricaud's modification,— that is 
described in the following subsection. 
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B. Paricaud's modified MIX1 theory (mMIXl) 

In the modification of the MIX1 theory introduced re- 
cently by Paricaud, 18 which we will refer to as mMIXl 
and ascribe a superscript mM, one keeps Eq. (|3.1[) , and 
hence Eqs. (|3.3l) - (|3.7p . except that Y^f is replaced by 



so that 
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By requiring Eq. (13. 12)) to be exact to first order in den- 
sity (third virial coefficient), i.e., cj^.y = Cfc ; y, one must 
have 



Yi 



mM 
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With this change Y?f -> Y^ M , Eq. dBHJ becomes 



a t 3 9 l j \P) = 
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Using the above results, the compressibility factor and 
Helmholtz free energy per particle in Hamad's proposal 
are given, respectively, by 

Z (P) = 1 +E-7^T ^ ■ (3-18) 
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and 

01 = -i+5>i» 

i 

where Z pure (y) and aP" ro (?/) are the compressibility fac- 
tor and the excess Helmholtz free energy per particle, 
respectively, of the one-component d-dimensional hard- 
sphere fluid at the packing fraction y. From Eqs. (|2.6[) 
and (|3.16p it follows that the fourth virial coefficient in 
Hamad's approximation is 

64 b 2 



(3.20) 



More in general, Eq. p,18[) yields 



n-2 



n— 2 



(3.21) 



may easily derive the Helmholtz free energy per particle, 
which turns out to be 

« SHY (^) b 3 B* 2 -b 2 B* 3 



k B T 



= -1 + 7^ Xj In (xjpXf) - 
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Also, Eq. (|3.24l) implies that the nth virial coefficient is 
given by 
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(3.26) 

while for the composition-independent fourth virial coef- 
ficients one gets the following explicit expression, 

rjSHY _ v d(fa - fa) / dn , Jd , _d R 

jM ~ 4(b 3 - b 2 ) ^ 1 j j u k ts ijt 

,d o \ V d(fa - fa) Id A o 

+a e B ijk ) - Rjh ^ - [ai^Bki 
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+af4B jt + af a fB jk +af4B il 
+<jj<jfB lk + a^afBij) . 



(3.27) 



D. The SHY proposal 

In Ref. [H we proposed the following ansatz for the con- 
tact values of the radial distribution functions, 
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where 
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This choice guarantees that gf ;) HY (p) is exact to first 
order in density and thus this approximation retains the 
exact second and third virial coefficients. When Eqs. 
(|3~22|) and (|3~23f are inserted into Eq. (|2~Tj) one gets 
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(3.24) 



where we have called B* = B n /{vd{<J )) n ~ l ; note that 
— > b n in the one-component limit. In Eq. (|3.24p we 
have expressed Z Y (p) — 1 as a linear combination of 
77/(1 — 7?) and Z pure (7?) — 1, with coefficients such that 
the second and third virial coefficients of the mixture are 
exactly reproduced. From the approximation (I3.24[) . one 



E. A nonlinear MIX1 theory 

As a final theoretical proposal, in this subsection we 
introduce a new extension of the MIX1 theory. 

The SHY approximation, Eq. ([3722]) . is a "local" ap- 
proximation with respect to density in the sense that 
the nonadditive contact value is expressed in terms of a 
reference contact value (here that of the one-component 
system) evaluated at precisely the same density. Anal- 
ogously, both the original MIX1 approximation, Eq. 
(|3.ip , and Paricaud's modified version, Eq. p. lip , can be 
termed "linearly non-local" since the nonadditive contact 
value is furthermore expressed in terms of the first den- 
sity derivative of the additive contact value. In contrast, 
Hamad's approximation, Eq. (|3.12p . is "nonlinearly non- 
local" because the reference contact value (again that of 
the one-component system) is taken at a totally different 
scaled density. 

Our nonlinear MIX1 (nlMIXl) approximation, labeled 
with nlM, is inspired in both Eq. (|33T|) and Eq. (|37l2]l . 
It consists of assuming that 



where 



™m 

j^-nlM _ ^ 1 *J 



1 1 ymM ' 

' ij 
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(3.29) 



Expanding in powers of X^ M - 1, Eq. (|3~2l?l) can be for- 
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mally rewritten as 
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(3.30) 



Comparison with Eq. p. lip shows that g™ M (p) can be 
seen as a first-order approximation of gf lM (p)- Using 
Eq. (|3.28p . the equation of state and Helmholtz free en- 
ergy per particle corresponding to the nlMIXl theory are 
given, respectively, by 

Z" 1M (p) = 1 + b 2 v d pY,XiXj49i? d (pX$ M ), (3-31) 
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where 



Q$*(p) ee v d (a d ) / dp'gtf d {p'). (3.33) 
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Note that, since gfj M (p) and gf 1M (p) coincide to first 
order in density, both give the same (approximate) third 
virial coefficient, namely Eq. with Y^f -s- F™ M . 

However, the mMIXl and nlMIXl theories differ at the 
level of the fourth virial coefficient. In this case, instead 
of Eq. dSU]) [with Y$ F™ M ] we have 
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It would be tempting to determine Xfj M in Eq. (|3T28|) 
by requiring agreement with the exact result to first order 
in density. This would give 



X)fc x kCk-ij 



(3.35) 



Unfortunately, however, this implies a wrong composition 
dependence of the higher order terms in the expansion of 
gij(p) in powers of p. In particular, 



XkXtCkl- 

k,e 



:ij 



K 2^k X k C k:ij / k<l 



While the left-hand side is quadratic in the mole frac- 
tions, the right-hand side is the ratio between a quartic 
function and a quadratic function. In order to avoid in- 
consistencies as in ()3.36j) we need Xf™ to be indepen- 
dent of the mole fractions. Apart from that, Xff can be 



freely chosen but we will keep the choice (|3.29[) in order 
to make contact with the mMIXl theory. 

Before closing this section, it is worth noting that, 
by construction, the nlMIXl theory is a priori not ex- 
pected to be accurate for strong negative nonadditivi- 
ties. This is because, on physical grounds, the param- 
eter Xf™ defined by Eq. ([3^|) must be positive def- 
inite. This in turn implies, from Eq. (|3.9p . the con- 
dition Aij > — (l — 2~ 1//d ). In the three-dimensional 
case, the above condition becomes Ajj > —0.21. As a 
matter of fact, the expansion (|3.30p does not converge 
if < — (l — 2~ 1 / d ). Notwithstanding this, from a 
practical point of view the nlMIXl theory keeps provid- 
ing meaningful results even if < — (l — 2~ 1 / d ), as 
will be seen in the next section. 



IV. RESULTS 

Thus far the development has been rather general in 
the sense that all the approximations we have discussed 
apply for any number of components N in the mixture 
and any dimensionality d. However, it is only formal un- 
less one specifies Z add (p), ag dd (p), and gf dd (p) in the case 
of all the MIX1 theories, and g puic (y), Z pure (y), a^ rc (y), 
and Ck-ij in the cases of Hamad's and the SHY approx- 
imations. In Ref. y we introduced for general d the fol- 
lowing approximation 



Ck\ij 



k;ij 
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0~ik 



Ti;jk&j;ik d- 
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(4.1) 



(4.2) 



This is exact when d = 1 and d = 3 and proved to be 
accurate also for d = 2. We will also use it here. 

As for the other remaining quantities, since the new 
numerical datai^ have been obtained for d = 3, we 
will restrict ourselves in the subsequent analysis only to 
this dimensionality. Therefore in the MIX1 theories we 
will take for Z &dd (p) and a% dd {p) the expressions given 
by the popular Boubhk-Mansoori-Carnahan-Starling- 
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Leland (BMCSL) equation of state ; 19 i 20 namely 
1 377 (a) (a 2 ) 
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while for gf dd (p) the choice will be the Boubhk-Grundkc- 
Henderson-Lee-Levesque (BGHLL) value o 19 ' 21 ' 22 given 
by 
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It follows from Eq. (|4.5[) that cj^ and c|^.^ are given by 
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Equation (|4.6|) is exact and agrees with Eq. (|4.1j) in the 
three-dimensional additive limit (63/62 = |, (Jkyij — > 0fc)- 
On the other hand, Eq. (|4.7p is approximate. According 
to Eq. (US), the quantity defined by Eq. ([3~33]) is given 
by 



ln(l-ry) + 3 
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Finally, in the case of the pure system, we will consider 
the expressions corresponding to the Carnahan-Starling 
(CS) equation of stated namely 



Z purc (y) 



1 + y + y 2 - y 3 
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With the above choices, the five approximations reduce 
to the CS equation of state in the one-component case 
<Tj = o\ In the additive limit, however, there are three 
independent proposals: BMCSL (to which the original 
MIX1 theories and its two variants, mMIXl and nlMIXl, 
reduce), Hamad's, and what we referred to as eCS in Ref. 
4. Of course, when nonadditivity is introduced, the five 
approximations differ from each other. 



A. Virial coefficients 

Figures Q~HB] show the comparison of the values of 
the composition-independent fourth virial coefficients, as 
given by the five theoretical proposals considered in this 
paper, with the recent data of Pellicane et aL 17 ^ 24 

One can immediately see that in the cases of -B1112 
and -B1222 the best overall performance is the one of the 
nlMIXl theory, followed closely by Hamad's approxima- 
tion. Also worth noting is that the mMIXl theory al- 
ready does a very good job, especially for the smaller size 
ratios, while the original MIX1 theory gives the poorest 
agreement. As far as -B1122 is concerned, the agreement 
of the theoretical predictions with the Monte Carlo data 
is much less satisfactory, getting poorer as the nonaddi- 
tivity parameter is increased. Here, no approximation 
is able to capture the negative values obtained by the 
Monte Carlo method for A > 0.2 and Hamad's approx- 
imation totally fails for small size ratios, irrespective of 
the value of the nonadditivity parameter. This is due 
to the fact that, while the four remaining theories cor- 
rectly reproduce the scaling behavior -B1122 ~ 0i02 m the 
high-disparity limit 02/01 0, Hamad's proposal yields 
B1122 ~ of in that limit. If one had to make a choice 
for this coefficient B1122, either the SHY proposal or the 
original MIX1 theory would perhaps be the ones to go 
for (especially for < A < 0.2 and 0.3 < er 2 /0i < 1), 
but with all due reserves. 

One might reasonably wonder whether the use of more 
accurate expressions for the additive contact values gf dd 
might correct the inability of the theories examined in 
this paper to predict negative values of the virial coeffi- 
cient -B1122 for small size ratio 02/01 and large nonad- 
ditivity parameter A. However, a closer analysis shows 
that this is not the case. According to Eq. (|2.6p . -B1122 oc 
3 C22;ii + 02 c ii;22 + 4er 3 2 ci2;i2- Therefore, at least one 
of the nonadditive second-order coefficients C22 ; ii, cn ; 22, 
and Ci2 ; i2 must be negative if B1122 < 0. In contrast, the 
additive coefficients c^ dd . e are positive for any 02/01 an d, 
as a consequence, all the approximate theories considered 
here predict positive values of Cij-ke for A > 0, as can be 
seen from Eqs. ([53]) . (j3"TTj) . ([3TTS]) . (j3T2"2"j) . and 
B. Compressibility factor 

To complement the above information on the virial co- 
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FIG. 1: Plot of the composition-independent fourth virial co- 
efficients B1112, -B1122, and -B1222 versus the size ratio o^/eri 
for a nonadditivity parameter A = 0.05. The dotted lines 
correspond to the original MIXI theory, Eq. (|3.7[1 . the short- 
dash lines correspond to the mMIXI theory, Eq. (|3,7[l with 
Y™ ->■ Y-j M , the thin solid lines correspond to the nlMIXl 
theory, Eq. (|3.34p . the long-dash lines correspond to Hamad's 
proposal, Eq. (|3.20[1 . and the thick solid lines correspond to 
the SHY proposal, Eq. |pT27) l. The symbols are Monte Carlo 
data from Ref. fl7l . 
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FIG. 2: Same as in[T] but for A = 0.1. 



efficients, in Figs.EHHwe present the results of our calcu- 
lations of the compressibility factors of binary nonaddi- 
tive hard-sphere mixtures and a comparison with avail- 
able simulation data. 

Figure [7] displays the dependence of Z on the nonaddi- 
tivity parameter (both positive and negative) for a sym- 
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a/a, 

FIG. 3: Same as inCfl but for A = 0.2. 



metric binary mixture at 7/ = 7r/30 ~ 0.105 and two val- 
ues of the mole fraction, namely x\ — 0.1 and x\ = 0.5. 
In this case, both the SHY proposal and the nlMIXl the- 
ory provide the best agreement, but the mMIXl theory 
also does a very good job. Hamad's proposal performs 
better at negative nonadditivities than at positive ones. 
As for the MIX1 theory, being linear in A, only captures 




0.0 0.2 0.4 0.6 0.8 1.0 

a 2 /o, 



FIG. 4: Same as in[T] but for A = 0.3. 



the region of small |A|. 

The superiority of Hamad's theory for negative non- 
additivities is confirmed by Fig. [5J which corresponds to 
the case of an equimolar asymmetric binary mixture with 
size ratio 02/ '(T\ — -| and a packing fraction 77 = 0.5. Here 
Hamad's approximation clearly outperforms all the rest. 
As a matter of fact, it becomes exact in the extreme limit 
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A — > — 14 A noteworthy feature is that, in contrast with 
both the original MIX1 and the mMIXl theories, the 
nlMIXl theory at least captures correctly the qualitative 
behavior of the compressibility factor with the nonaddi- 
tivity parameter for negative values and, in particular, 
the initial decay. This is remarkable in view of the fact 
that, as discussed at the end of the preceding section, the 
nlMIXl theory is not expected to hold if A < -0.21. 



Finally, in Fig. [9] we present the results obtained for 
the size-ratio dependence of the compressibility factor 
for r\ = 0.2, a positive nonadditivity A = 0.2, and two 
compositions. In agreement with the behavior observed 
in Fig. [7] for A > 0, we see from Fig. U that the SHY is 
the superior theory also in the asymmetric case, although 
all the theories, with the exception of the MIX1, tend to 
coincide as the asymmetry increases. It is noteworthy 



10 



N 



2.2 
2.0 
1.8 
1.6 
1.4 
1.2t 



MIX1 


x=0.5// ."' 


mMIXl 


i //• • 

y ' 


" nlMIXl 


/ •' /: 


H 


/ 


emir 

brl I 




• MC 









-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

A 

FIG. 7: Plot of the compressibility factor Z versus the non- 
additivity parameter A for a symmetric binary mixture of 
nonadditive hard spheres at r\ = n/30 and two different 
compositions. The dotted lines correspond to the original 
MIX1 theory, Eq. (|3.3[) . the short-dash lines correspond to 

Y^ M , the thin 



the mMIXl theory, Eq. I03J with Y^f 



solid lines correspond to the nlMIXl theory, Eq. (|3.31[) . the 
long-dash lines correspond to Hamad's proposal, Eq. (|3.18[) . 
and the thick solid lines correspond to the SHY proposal, Eq. 
(|3.24[) . The symbols are results from Monte Carlo simulations 
(Refs.QMl]). 
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FIG. 8: Plot of the compressibility factor Z versus the non- 
additivity parameter A for an equimolar asymmetric binary 
mixture of nonadditive hard spheres with size ratio 02/ &i = § 
at 77 = 0.5. The dotted line corresponds to the original 
MIX1 theory, Eq. (|3.3[) . the short-dash line corresponds to 

Y^ M , the thin 



the mMIXl theory, Eq. $3^ with Y™ 



solid line corresponds to the nlMIXl theory, Eq. H3.31|) . the 
long-dash line corresponds to Hamad's proposal, Eq. I|3.18[l . 
and the thick solid line corresponds to the SHY proposal, Eq. 
(|3.24p The symbols are results from Monte Carlo simulations 
(Refill). 



that both the mMIXl and the nlMIXl theories do a very 
reasonable job, better than Hamad's proposal. 

C. Demixing 

The availability of analytical expressions for the 
Helmholtz free energy per particle a in all the previ- 
ous theories [cf. Eqs. (03), (|3T9| . ([3^5]) . and ([3732]) ] 
may be exploited to address the problem of demixing in 
mixtures with positive nonadditivity. For simplicity, we 
will restrict ourselves here to binary mixtures. Since in 
these systems the temperature only plays the role of a 
scaling factor and a spinodal instability occurs, the mix- 
ture will phase separate into two liquid phases (I and II) 
of different composition x\ and ij 1 . For given size ra- 
tio o^joi and nonadditivity A, by equating the pressure 
(p 1 = p n ) and the two chemical potentials (p\ = /i* 1 , 
p\ = p\ l ) of both phases, one may obtain p 1 , p u , and a:* 1 
as functions of x\ and thus derive the coexistence curve 
in the p-x\ plane. The chemical potentials are defined by 
Pi = d(pa)/dpi. In the binary case, this is equivalent to 



P P 
p 1 = a-\ \-(l-xx)a x , p 2 =a-\ xia x , (4.12) 

P P 



where a x = (da/dxi) p . The two branches I and II of 
the coexistence line meet at the critical consolute point 



(p c , x\ c ), which can be determined by the two conditions 



(4.13) 



= ( a pp + -a p I a xx - a xp , 



3a x 



3 1 a 

G 



tj Xp 



®>xp 



(4.14) 



Here, as in Eq. (14.12[) . each subscript x or p represents 
a derivative with respect to x\ or p, respectively. For 
symmetric mixtures, the critical composition is fixed, 
xic = 0.5. 

In Fig. [10] we display the behavior of the reduced crit- 
ical density p* — Pc^s m symmetric mixtures, where 
°eff = Si j x i x j a iji as a function of the nonadditivity 
parameter A for Hamad's theory, the SHY proposal, the 
mMIXl and nlMIXl theories, and the available simula- 
tion data. The original MIX1 theory has not been in- 
cluded since it has already been proved that it yields a 
poorer performance than Hamad's theory which is the 
least accurate in this instance. Note that all theoretical 
results underestimate p* and are very close to one an- 
other with perhaps a slightly better overall performance 
of the mMIXl and the SHY. The use of the effective 
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FIG. 9: Plot of the compressibility factor Z versus the size 
ratio oil(T\ for binary mixtures of nonadditive hard spheres 
with A = 0.2 and x\ = 0.75 (upper panel) and x\ = 0.5 
(lower panel). The dotted lines correspond to the original 
MIX1 theory, Eq. (|3.3(l . the short-dash lines correspond to 

Y^ M , the thin 



the mMIXl theory, Eq. ((OJ with Y™ 



solid lines correspond to the nlMIXl theory, Eq. (I3.31[) . the 
long-dash lines correspond to Hamad's proposal, Eq. (13.181) . 
and the thick solid lines correspond to the SHY proposal, Eq. 
(|3.24[) The symbols are results from Monte Carlo simulations 
(Refill). 



diameter a e s to define the reduced critical density in 
Fig. [TU] is motivated by the fact that p* is well defined 
for high nonadditivities, including the Widom-Rowlinson 
limit (<7i — i72 *C cri2 or A — >• oo). 

As far as the liquid-liquid coexistence curve is con- 
cerned, this may be represented in different thermody- 
namic planes. Here we have chosen the pa\-x\ and the 
pa^/ksT-xi planes. Further, given the previous analy- 
sis concerning the comparison of the theoretical critical 
consolute points and simulation results, and the techni- 
cal difficulties associated with the actual computation of 
the coexistence curves, only the results for the SHY, the 
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FIG. 10: Plot of the reduced critical density p* — pc&ca versus 
the nonadditivity parameter A for symmetric binary mixtures 
of nonadditive hard spheres. The short-dash line corresponds 
to the mMIXl theory, the thin solid line corresponds to the 
nlMIXl theory, the long-dash line corresponds to Hamad's 
proposal, and the thick solid line corresponds to the SHY pro- 
posal. Th e sy mbols are results from Monte Carlo simulations 
(Refs.|25H27|). 



mMIXl and nlMIXl theories will be presented. A com- 
parison of available simulation results for liquid-liquid 
coexistence is done both for symmetric and asymmetric 
mixtures in Fig. 111! where the theoretical critical con- 
solute points have also been included. Notice that the 
qualitative trends observed in the simulations are well 
captured by all the theoretical developments, but in all 
instances they tend to underestimate the actual values of 
the reduced pressure and the reduced density along the 
coexistence. In particular, all theories correctly predict 
that the demixing transition occurs for lower densities 
as the nonadditivity parameter increases. Moreover, at 
a fixed value of A the coexistence densities (if measured 
in units of the diameter of the smaller spheres) decrease 
with increasing size asymmetry. Similar trends are ob- 
served for the pressure. On the quantitative side, par- 
ticularly in the density vs composition plane, albeit not 
very accurate, the SHY outperforms the other theoretical 
approximations. 

V. CONCLUDING REMARKS 

In this paper we have provided a self-contained presen- 
tation of different theoretical developments to describe 
the thermodynamic properties of nonadditive hard-core 
mixtures. In particular, complementing the effort ini- 
tiated in our previous paper on this subject,— apart 
from repeating the SHY proposal and the extension of 
Hamad's approach to general dimensionalities, here we 
have provided extensions of the original MIX1 and Par- 
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FIG. 11: Liquid-liquid coexistence curves for several binary 
mixtures of nonadditive hard spheres in the reduced density 
peri versus composition xi plane (top panel) and in the re- 
duced pressure pa^/ksT versus composition Xi plane (bot- 
tom panel). From top to bottom, the set of curves corre- 
spond to (a"2/o"i, A) = (1,0.1) (absent in the bottom panel), 
(1,0.2), and (5/6,0.1818). The short-dash lines correspond 
to the mMIXl theory, the thin solid lines correspond to the 
nlMIXl theory, and the thick solid lines correspond to the 
SHY proposal. The diamonds indicate the locations of the 
respective critical consolute points. The other symbols are 
results from Monte Carlo simulations: Ref. [25| (filled circles), 
Ref. [13 (open circles), and Ref. @ (filled squares). 



icaud's modified MIX1 (mMIXl) theories valid for all d. 
We have introduced as well a new nonlinear extension 
of the MIX1 (nlMIXl) theory, also valid for arbitrary 
d. In all instances, explicit expressions have been pro- 
vided for the contact values of the radial distribution 
functions, the compressibility factor, the Hclmholtz free 
energy, and the second, third, and fourth virial coeffi- 
cients. The expressions for gij(p) and Z(p) are given in 
terms of either gf dd (p) and Z a,dd (p) in the case of all the 
MIX1 theories, or in terms of g puTe (y), or equivalently of 
ZP UTC (y) = l+2 d - 1 ?; ff P uro (y), in the cases of Hamad's and 



the SHY approximations. For the sake of illustration and 
restricting to three-dimensional systems (d = 3), we have 
taken as input the BMCSL equation of state for Z add (p) 
and the BGHLL contact values for gf dd (p) in the MIX1 
theories, and the CS equation of state for Z puYe (y) in the 
SHY and Hamad's proposals. 

To our knowledge, the idea of starting from the con- 
tact values of the radial distribution functions in the case 
of the MIX1 theories has not been considered before. 
This allowed us to construct the nonlinear extension. Of 
course, while in the case of mixtures the compressibil- 
ity factor is determined uniquely once the contact values 
of the radial distribution function are given, the recipro- 
cal is not true. Hence, the expressions we have provided 
for these contact values are a further contribution of this 
work. 

We have carried out three kinds of comparison between 
the five theories and "exact" numerical results. First, the 
theoretical predictions of the composition-independent 
fourth virial coefficients have been tested against new 
available Monte Carlo dataJ^ In the cases of -B1112 and 
-B1222, the best overall agreement with the Monte Carlo 
values are obtained with the nlMIXl theory, followed by 
Hamad's proposal. As for -B1122, none of the theories does 
well at high asymmetry and nonadditivity, the discrep- 
ancies being especially important in the case of Hamad's 
approximation. 

As is well known, the first few virial coefficients are rel- 
evant to the equation of state in the low-density regime 
but not generally beyond it. Thus, in order to test the 
theoretical approaches at finite densities, we have made 
use of available simulation data for the compressibil- 
ity factor i 10 ' 11 ' 13 The emerging scenario is that Hamad's 
approximation is excellent for negative nonadditivities, 
while the SHY proposal is the preferable one for positive 
nonadditivities. 

Within the limited set of compressibility factors that 
we have analyzed, it is fair to say that the new nlMIXl 
theory proposed in this paper is rather satisfactory and 
seems to be a good compromise between accuracy and 
simplicity. Further assessment of this assertion is pre- 
cluded at this stage due to the scarcity of the data. One 
of our hopes is therefore that the present paper may en- 
courage more work on the subject. 

Finally, the critical behavior and liquid-liquid coexis- 
tence of nonadditive hard-sphere mixtures with positive 
nonadditivity has been examined. While the quantita- 
tive agreement is not satisfactory, all theories seem to 
capture correctly the qualitative trends obtained in the 
simulation. In this case our original SHY proposal gives 
the best performance, but again the limited availability 
of data prevents us from carrying out a more thorough 
analysis. Once more we hope that our findings may lead 
to the further needed work on this matter. 
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